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The neutron is largely spherical and incompressible in atomic nuclei. These two properties are 
however challenged in the extreme pressure environment of a neutron star. Our variational compu- 
tation within the Cornell model of Coulomb gauge QCD shows that the neutron (and also the A 3 ' 2 
baryon) can adopt cubic symmetry at an energy cost of about 150 MeV. 

Balancing this with the free energy gained by tighter neutron packing, we expose the possible soft- 
ening of the equation of state of neutron matter. 



The atomic nucleus is a close-packed system analogous 
to a cold quantum liquid. Its hadronic components are 
very incompressible, as revealed by the nuclear radius 
scaling with the mass number as ta — r^A^ 1 / 3 , with 
unit of length r$ ~ 1.2 fm. Nucleons there are quite op- 
timally packed, as the (charge) nucleon radius has been 
precisely determined [l| in electron scattering, and is ap- 
proximately 0.88 fm. If we (grossly) think of the nucleon 
as a solid sphere, 74% of the nuclear volume is occupied, 
near the Keplerian limit Q reached by hep or fee lattices. 

In neutron stars, gravity-compressed nuclear matter 
finds extreme pressure environments that changes some 
of its properties. Recently Q a two-solar mass neutron 
star has been discovered that takes nuclear physics to 
the limit (and we have recently proposed [4| that this is 
so much so as to constrain gravity for strong field inten- 
sity). At high pressure and density in the core of the 
star many exotic phases of hadron matter have been pro- 
posed, including a neutron superconductor/superfluid, 
various meson condensates, quarkyonic or Color-Flavor- 
Locked condensates, etc. Some of them have already 
been ruled out by the new, superheavy neutron star 
At the crust one expects a thin atomic sheet @ followed 
by neutron-rich nuclear matter, and finally almost pure 
neutron matter (with a small number of protons in /3- 
cquilibrium) [jj. 

One intermediate phase that has long been pro- 
posed [H is a neutron crystal. In analogy with condensed 
matter systems (saliently 3 He, a quantum liquid through 
T = that solidifies upon compression) neutron matter is 
assumed to locally adopt a periodic lattice-like arrange- 
ment to maximize density. 

Our new observation is that when the environment looses 
its isotropy, the nucleon (largely spherical at rest) is sub- 
ject to directional stresses and may deform to minimize 
the evacuated volume. Thus, for mass-energy densities e 
larger than 140 MeV j fm 3 ~ 3to4 neutrons can no more 
be considered pointlike objects subject to the rules of 
local, rotationally-invariant quantum field theory as ele- 
mentary fields. Their composition and structure should 
start being taken into account, perhaps revisiting some of 



the hypothesis underlying the nuclear many body prob- 
lem approach based on QFT methods (particularly, ro- 
tational invariance). 

We report an estimate of the deformability of the neu- 
tron with the renowned Cornell Hamiltonian, 

H{x) = &(x){-ia- V + to/3)* + (1) 

d 3 y^(x)T a y(x)V(x-y)^(y)T a y(y) 

a field theory model of Quantum Chromodynamics, 
tractable by many-body methods common in nuclear and 
condensed matter theory. The relativistic quark fields 
* interact via the color charge density ^(x)T ai i l (x). 
The model has been am ply exploited to treat the heavy 
quarkonium spectrum ^qL 1 1 Ol] , as well as both light mesons 
and baryons HJ El with the Cornell potential V(r) = 
or — 

The model ground state ( Po quark-antiquark con- 
densed vacuum) is treated in BCS approximation to gen- 
erate constituent quarks from the current to — 5 MeV 
quarks in the Hamiltonian, then the second-quantized 
wavefunction appropriate for 3-quark baryons in terms 
of Bogoliubov-rotated quark creation operators is 

\n) = — p^ iA2A3 (fci,fc 2 ,fc 3 )P±" Bi b+ in> . 

17 6 N feiAi* k 2 \ 2 j k 3 \ 3 k< 

(2) 

To approach lowest-lying baryons by the Rayleigh-Ritz 
variational principle, we employ the separable wavefunc- 
tion (p, A being 3-body Jacobi coordinates with Cartesian 
components p x , X z , etc.) 
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To choose a convenient variational wavefunction we ob- 
serve that the graph of (k% + ky + k^)~ = 1, which 
is depicted in figure [T] for even N, interpolates between 
spherical symmetry for N — 2 and octahedral (the 
cube's) symmetry, reached as N — > oo. Therefore one 
can take an appropriately normalized ipN from a pre- 
cise numerical solution of the two-body problem and use 



FIG. 1. Trial wavefunction that interpolates between sphere (for N = 2), and cube (as TV — > oo) for N = 2,4, 8, 12. 



ipN (fax + ^ + to study the three-quark sys- 

tem transiting between both symmetry groups. 

The Hamiltonian expectation value (N\H\N) is a func- 
tion of the two variational parameters a p , a\, and is min- 
imized respect to them to find the best approximation 
to the neutron mass within this family of functions for 
N = 2. To compute the matrix elements we perform all 
relativistic spinor sums numerically, and the three- body 
nine-dimensional integral by Monte Carlo methods [13j . 
The model approximation to the nucleon mass is 980 
MeV (with a Monte Carlo error of 40 MeV), the phys- 
ical mass being 940 MeV . Should higher precision be 
necessary one should resort to lattice gauge theory [l4j . 
This is not so for our purposes and we proceed to change 
the wavefunction symmetry, repeating the minimization 
for varying N . 

We find that the neutron mass increases by about 
150 MeV between N = 2 and N = 18 (the last already 
very close to a cubic neutron) , as depicted in figure [2 

As a check of the spin-independence of our result, we 
have repeated the calculation for the A(1232) baryon, 
where all three quark spins are parallel, and obtained an 
equal excitation energy. The similar outcome is reported 
in figure [3] 

We can offer the following heuristic argument sup- 
porting these numerical findings. If one takes a high 
N (close to cubic) variational wavefunction, and per- 
forms a spherical harmonic analysis, 75% of the prob- 
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• Deformed neutron 

□ Experimental P=+ resonances 
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FIG. 2. Open squares: experimental neutron mass and known 
positive parity resonance excitations [Tq| . Filled circles: vari- 
ational minimum for the nucleon mass with the wavefunction 
in Eq. @. At a cost of 50 MeV the neutron already looses 
spherical symmetry, adopting an almost totally cubic shape 
for 150 MeV. 



ability resides in the loo spherically symmetric state, 
and an additional 24% in one of the d-wave Y2m (mostly 
to = 0) states. Weighing the experimental masses (sev- 
eral positive-parity excitations of the nucleon are also 



shown in figure [2]) of ground J 1 



i + 



(940 MeV) and 



d-wave excited J p = § (1680 MeV) nucleons in the 
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FIG. 3. Open squares: experimental A-baryon mass and 
known positive parity resonance excitations 15j. Filled cir- 
cles: variational minimum for the A mass with the wavefunc- 
tion in Eq. (|3}. The excitation cost to abandon spherical 
symmetry is very similar to that reported in figure [5] for the 
neutron. 
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FIG. 4. Mass of nucleon and A-baryon as function of only 
one variational parameter ct p = ct\. In addition to a mini- 
mum corresponding to the best Rayleigh-Ritz representation 
of the particles within the function family, the energy cost of 
compression or dilatation of the wave function (at constant 
shape) can be appreciated. 



proportion 3:1 yields 1125 MeV. This entails a cost of 
about 180 MeV to change the neutron symmetry from 
spherical to cubic. Thus, from a knowledge of the experi- 
mental spectrum and symmetry considerations alone one 
can already obtain the approximate answer. The same 
reasoning applies to the 150 MeV excitation energy of 
the A, that can likewise be interpreted since the relevant 
d-wave resonance energy is tagged by the A 5 / 2+ (1900) 
and the A 7 / 2 +(1950). 

Although deformability at constant hadron volume, 
and not compressibility of the hadron, was the goal of 
our investigation, we can provide a feeling for the latter 
by setting both Jacobi scales equal a p — ot\ and studying 
the mass as function as this variational parameter alone. 
We refer the reader to figure |4l A variational minimum 
is clearly seen for both nucleon and A, corresponding 
to their best Rayleigh-Ritz representation. Compress- 
ing or dilating the neutron can be faked by varying this 
scale, and the cost of doing so is evident in the figure, 
and of the order of magnitude that one naively expects. 
The neutron seems to be somewhat stiffer to momentum- 
space compression (position-space dilatation), while both 
particles react similarly to momentum-space dilatation 
(position-space compression). 

The excitation energy that we have found is almost 
equal to that necessary to produce hyperons, since the 
lightest strange particle, A, has a mass of about 1115 
MeV . Thus, if the nucleon's Fermi momentum is high 
enough to relax some of the d quarks to s quarks, it is also 
high enough to deform the nucleon, and both phenom- 
ena appear simultaneously. Their interplay is beyond our 
current undertaking. 

At high enough pressure, the 150 MeV energetic cost 
to deform the neutron can be provided by the Helmholtz 
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FIG. 5. Interstitial volume gain AV by deformation of the 
neutron wavefunction at constant neutron volume for two 
cubic lattices. The gain is naturally smaller for the more 
packing-efficient face-centered-cubic lattice, but still signifi- 
cant. 

free energy PAV obtained upon reducing the interstitial 
volume between neutrons thanks to the cubic shape. The 
maximum gain as N — > oo in an fee lattice, with pack- 
ing fraction 0.74, is 26% of the initial volume occupied 
by a spherical neutron. Estimating this volume for ar- 
bitrary N through simple three-dimensional integration 
(see figure [5]) we can, at fixed pressure, obtain the conse- 
quent increase in density. We plot in figure [6] the effect 
on the equation of state of neutron matter of the neu- 
tron deformability that we have found. The equation of 
state becomes softer (smaller pressure at given density) 
at densities relevant for neutron star structure. 

To conclude, let us observe that the bulk proper- 
ties of neutron matter in the star will be affected by 
the neutrons adopting cubic symmetry (and patches of 
neutron matter becoming crystalline), particularly the 
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FIG. 6. Effect on the state of the art equation of state[lq.ll7l] 
of neutron matter of the deformability of the neutron. The 
entire interval shown is relevant for neutron star structure. 
Nucleon deformability produces a softening of the equation 
of state additional to other known sources. 



Young modulus, the moment of inertia, etc. and this will 
surely provide detailed tests for future work. 
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